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Abstract 



We analyze daily prices of 29 commodities and 2449 stocks, each over a 

period of ~ 15 years. We find that the price fluctuations for commodities 

have a significantly broader multifractal spectrum than for stocks. We also 

propose that multifractal properties of both stocks and commodities can be 

attributed mainly to the broad probability distribution of price fluctuations 

and secondarily to their temporal organization. Furthermore, we propose 

that, for commodities, stronger higher order correlations in price fluctuations 

result in broader multifractal spectra. 
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The study of economic markets has recently become an area of active research for physi- 
cists [1], in part because of the large amount of data that can be accessed for statistical 
analysis. Markets are complex systems for which the variables characterizing the state of 
the system — e.g., the price of the goods, the number of trades, and the number of agents, 
are easily quantified. These variables serve as fundamental examples of scale invariant 
behavior — the scaling laws are valid for time scales from seconds to decades. 

Much interest has concentrated on stocks, where a number of empirical findings have been 
established, such as [2] (i) the distribution of price changes is approximately symmetric and 
decays with power law tails with an exponent a + 1 fti 4 for the probability density function; 
(ii) the price changes are exponentially (short-range) correlated while the absolute values of 
price changes ( "volatihty" ) are power-law (long-range) correlated. 

Unlike stock and foreign exchange markets, commodity markets have received little at- 
tention. Recently, it was found [3] that commodity markets have qualitative features similar 
to those of the stock market. This similarity is intriguing because the commodity market 
has special features such as: (i) most commodities require storage; (ii) most commodities 
require transportation to bring them to the market from where they are produced; and (iii) 
it is plausible that commodities may exhibit a slower response to change in demand because 
the price depends on the supply of the actual object. 

The multifractal (MF) spectrum [4] reflects the n-point correlations [5] and thus provides 
more information about the temporal organization of price fluctuations than 2-point corre- 
lations. Previous work reports a broad MF spectrum of stock indices and foreign exchange 
markets [6-12]. Two recent models [13,14] explain the observed MF properties by assum- 
ing that price changes are the product of two stochastic variables, one being uncorrelated 
and normally distributed and the other being correlated and log-normally distributed. The 
price changes predicted by these models do not have the power law probability distribution 
[2, 3] observed empirically, and thus shuffling the price changes signiflcantly narrows the MF 
spectrum. 

Here we show that the MF properties of commodities and stocks result partly from the 
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temporal organization and partly from the power law distribution of price fluctuations. We 
also conjecture that it is feature (iii) of commodity markets that leads to a broader MF 
spectrum of commodities compared to stocks. 

We define the normalized price fluctuation ( "return" ) as 

where here At — 1 day, S{t) is the price, and a is the standard deviation of ln5'(i + At) — 
In S'(t) over the duration of the time series (typically 15 years). We use the multifractal 
detrended fluctuation analysis (MF-DFA) method [15] to study the MF properties of the 
returns for stocks and commodities. Following [15], given a time series Xk we execute the 
following steps: (i) Calculate the profile; Yi = J2k=i[^k — {x)] i — 1,...,N, where N 
is the length of the time series, and (x) is the mean, (ii) Divide Yi into Ns = int{N/s) 
segments, (iii) Calculate the local trend yu{i) for segments u — l,...,Ns by least square 
polynomial fit. (iv) Determine the mean square fluctuation in each segment ^2(5, = 
-,J:t=l{\y{^-l}s+^ - y.m^- (v) Evaluate F,{s) = ^ E^=i ^^2(5, i^)^/'- The scaling function 
of moment q, Fq{s) [15] follows scahng law 

F,{s) s^^'^l (2) 

Negative q values weight more small fluctuations, while positive values of q give more weight 
to large fluctuations. 

When the contribution of the small fluctuations is comparable to the contribution of the 
large fluctuations, the series is monofractal and T{q) — hq where h — const is the Hurst 
exponent. If r(g) nonlinearly depends on q, the series is MF. The Legendre transform of 
T{q) is 

f(h) = qh-T(q) where h=^^. (3) 

aq 

Monofractal series have only one value of h, unlike MF series which have a distribution of h 
values. 



We analyze a database consisting of daily prices of 29 commodities [16] and 2449 
stocks [17] spread over time periods ranging from 10 to 30 years (the average period is 
a; 15 years). Figure 1 shows the price of a typical commodity, corn, and a typical stock, 
IBM, and their corresponding returns. One striking difference between the commodity and 
the stock is that the commodity returns are more "clustered" into patches of small and 
large fluctuations. This feature is not reflected in the distribution and the autocorrelation 
function of the price fluctuations [3] . 

To emphasize the clustering of commodities, we shuffle the returns by randomly ex- 
changing pairs, a procedure that preserves the distribution of the returns but destroys any 
temporal correlations. Specifically, the shuffling procedure consists of the following steps 

(i) Generate pairs (m, n) of random integer numbers (with m,n < N) where N is the 
total length of the time series to be shuffled. 

(ii) Interchange entries m and n. 

(iii) Repeat steps (i) and (ii) for 20A'" steps. (This step ensures that ordering of entries in 
the time series is fully shuffled.) 

To avoid systematic errors caused by the random number generators, the shuffling proce- 
dure is repeated with different random number seeds for each of the 2449 stocks and 29 
commodities. The shuffled commodity series of necessity loses its clustering [18] [Fig. 1(e)]. 
On the other hand, the shuffled stock series resembles the original one. 

First we compare the 2-point correlations using DFA [19] of the shuffled and the un- 
shuffled returns for commodities and stocks (Fig. 2). Both corn (and all other commodities 
analyzed) and IBM (and all other stocks analyzed) have returns uncorrelated for time scales 
larger than one day [2, 3]. Thus, studying the 2-point correlations is not sufflcient to uncover 
the clustering of the commodity returns. 

Next, we analyze the MF properties of the returns of stocks and commodities. Figure 3(a) 
displays separately the averages 
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1 ^ 

^av(g) = V^E^i(5)' (4) 
i=l 

for = 29 commodities, and for N = 2449 stocks. Note that (i) the scahng exponents, 
Tav(?)|g<o for commodities and stocks significantly differ, whereas Tav(?)|g>o are similar, 
suggesting that commodities are similar to stocks for the large fluctuations and they differ 
for the small fluctuations [20]; (ii) we observe that after shuffling the returns, Tav{q) for stocks 
hardly changes for g < 0, but Tav(g) for commodities changes and becomes comparable to 
stocks for the entire range of q [21]. 

In order to study the contribution of the power law tails of the returns on the MF 
spectrum we generate surrogate data sets (i) with a normal distribution and (ii) with power 
law tails with a 3 (as is observed empirically [2,3]). Figure 3(b) displays Tav(?) averaged 
over 2449 realizations of surrogate data, each with 3000 data points. The Tav{q) of the 
surrogate power law distributed data is very close to the Tg_y{q) of stocks after shuffling. 
This indicates that a significant part of the Tav(?) spectrum of stocks and commodities 
comes from the power law distribution of the returns. Note that there is a small difference 
in Tscvlq) of stocks before and after shuffling, indicating that the power-law distribution 
of the returns is not the only source of multifractality but that there is also a relatively 
small contribution due to temporal organization of returns. For commodities this temporal 
organization is more dominant. 

Figure 4 displays the MF spectrum of the unshuffled and shuffled returns for commodities 
and stocks. The temporal organization of commodity returns is reflected in the fact that 
the MF spectrum for unshuffled commodities is broader than for shuffled commodities and 
stocks. 

Demand fluctuations drive price fluctuations, and it is plausible that stocks respond more 
quickly than commodities to demand changes. Stochastic perturbations, together with the 
immediate price response to demand changes, may weaken the existing higher order temporal 
organization, which may be the reason for less clustering for stocks (Fig. 1). Commodities, 
on the other hand, have a slower response. Thus, small or short-time perturbations are 
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"felt" less by commodities than by stocks. This scenario is consistent with the appearance of 
patches of small commodity fluctuations followed by patches of large commodity fluctuations 
(Fig. 1). We conjecture that the more homogeneous returns of stocks explain the difference 
between the MF properties of stocks and commodities. 

In summary, we find that commodities have a broader MF spectrum than stocks. A 
major contribution to multifractality is the power law tail of the probability distribution 
of the returns. Moreover, the MF spectra of stocks and commodities are partly related to 
the power law probability distribution of returns and partly to the higher order temporal 
correlations present. 
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FIG. 1. Analysis of 2048 daily returns covering the time period May 1993 - June 2001. Price 
of (a) corn, a typical commodity, and (b) IBM, a typical stock. Returns for (c) corn and (d) 
IBM; unlike stocks, commodity returns appear more clustered into "patches" of small and large 
fluctuations, even though there are no 2-point correlations (Fig. 2). Returns for (e) corn and 
(f) IBM after shuffling the returns; because all n-point correlations are now removed, stocks and 
commodities look similar and do not appear to cluster into "patches." 



10 



cj Typical Commodity, Com Typical Stock, IBM 



§ 10° 
2 



.Mll| 1 




1 ' ' 

•=1 


■ (a) 










0.5 




• returns 


\ * 


□ shuffled returns 


10' 


10^ 


10^ 



- (b) 








" H 




10' 


10" 10' 



s (days) 



FIG. 2. DFA analysis of two-point correlations of returns of (a) a typical commodity (corn) 
and (b) a typical stock (IBM), before and after shuffling. The exponent of 0.5 indicates that both 
the commodity and the stock are uncorrelated in time, so the two-point correlation function does 
not provide information regarding the clustering into patches (Fig. 1). 
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FIG. 3. (a) Ta,v{q) for returns and shuffled returns for 29 commodities and 2449 stocks. To 
better visualize the results we plot rav(q') — q'/2 instead of Ta,y{q). The exponents Ta,v{q) are calcu- 
lated for window scales of 10 — 100 days. After shuffling T^viq) are comparable for both stocks and 
commodities, (b) Ta,y{q) spectrum of the returns and shuffled returns for stocks, compared with 
uncorrelated surrogate data with Gaussian probability distributions and power law probability dis- 
tributions (with power law exponent a « 3). After shuffling, Ta.y{q) for stocks becomes comparable 
with Ta,y{q) of the surrogate data obtained for the power law probability distribution. 
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FIG. 4. fa.v{h) vs. h computed from the Ta,y{q) vs. q for commodities and stocks. Shown are also 
both unshuffled and shuffled returns. For commodities, the unshuffled returns show a broader MF 
spectrum than for the shuffled returns, consistent with the hypothesis that the broad MF spectrum 
arises due to temporal organization and the broad power law distribution of price fluctuations (cf. 
Fig. 3b). 
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